Laser induced breakdown of the magnetic field reversal symmetry in the propagation 

of unpolarized light 
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We show how a medium, under the influence of a coherent control field which is resonant or 
close to resonance to an appropriate atomic transition, can lead to very strong asymmetries in the 
propagation of unpolarized light when the direction of the magnetic field is reversed. We show how 
EIT can be used to mimic effects occurring in natural systems and that E1T can produce very large 
asymmetries as we use electric dipole allowed transitions. Using density matrix calculations we 
present results for the breakdown of the magnetic field reversal symmetry for two different atomic 
configurations. 

PACS No(s).:42.50.Gy,78.20.Ls,33.55.-b 

I. INTRODUCTION 

It is well known how an isotropic medium becomes 
anisotropic by the application of a magnetic field |Q. In 
the special case when the magnetic field is parallel to the 
direction of the applied field and if we include the elec- 
tric dipole contribution to susceptibilities, then the two 
circularly polarized light waves travel independently of 
each other. The propagation itself is determined by the 
magnetic field dependent optical susceptibilities \±- If 
we write the incident field of frequency to in the form 



h = -[\e 2 * m ^\ 2 + \e 27!lkl *-\ 2 } 

= T - [exp{-47rfcZIm(x+)} + exp{-47rfcZIm(x_)}] . (6) 

Thus the output intensity is a symmetric function of x+ 
and X-- If the susceptibilities \± obey the following re- 
lation when the direction of the magnetic field is reversed 



X±(B) = X T (-B) 



then 



E= (£ + e+ + £_e_)e 



ikz — iut 



c.c. 



(1) 



I (B)=I (-B) 



(7) 



(8) 



where 



x ± iy 



£± 



V2 J V V2 

the output field is given by 



i/o = ooe + c.c, 



where 



So = £ + e+e 2 ™ klx+ 



-E-€-e 2mklx - 



(3) 



(4) 



The susceptibilities x± a l so depend on the frequency of 
the applied field. The rotation of the plane of polariza- 
tion and the dichroism can be calculated in terms 
of the real and imaginary parts of \± ■ An interesting sit- 
uation arises if the incident pulse is unpolarized. In that 
case there is a random phase difference between £ x and 
£ y and the intensities along two orthogonal directions are 
equal, i.e., 



(£*£ x ) = (£* y £ y 



I 

2' 



(£*£ y ) = 0, 



(5) 



/ being the intensity of the incident pulse. From Eqs. 
we can evaluate the output intensity I = (|£o| 2 ): 



Thus for unpolarized light the output intensity is the 
same whether the magnetic field is parallel or antiparal- 
lel to the direction of propagation of the electromagnetic 
field as long as (0) is satisfied. 

In this paper we investigate if the transmission of unpo- 
larised light through an otherwise isotropic medium can 
be sensitive to the direction of the magnetic field. We 
demonstrate how a suitably applied control field could 
make the transmission dependent on the direction of the 
magnetic field. This is perhaps the first demonstration of 
the dependence of transmission of the unpolariscd light 
on the direction of B. For a large range of parameters 
we find that transmission could be changed by a factor of 
order two. We also report a parameter domain where the 
medium becomes opaque for one direction but becomes 
transparent for the reversed direction of the magnetic 
field. This work is motivated by the phenomena of optical 
activity and the magneto-chiral anisotropy which occur 
in many systems in nature fg-[l(|. The latter effect has 
recently become quite important. Several ingenious mea- 
surements of this effect have been made though the effect 
in natural systems is quite small. The smallness of the 
effect arises from the fact that the effect involves a combi- 
nation of electric dipole, magnetic dipole and quadrupole 
effects [11- 14]. The magnetochiral anisotropy is just the 
statement 1(B) not equal to I(—B). The effect wc report 
is analogous but quite different in its physical content. 
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In this paper we show how rather large asymmetry 
between Iq(B) and Iq{—B) can be produced by using a 
coherent control field. The asymmetry could be large as 
we use only the electric dipole transitions. To demon- 
strate the idea we consider different specific situations 
depending on the transition on which the control field is 
applied. The control field can be used to modify, say, 
X+ leaving X- unchanged |p|,^|j6|JT5|l . Thus in presence 
of the control field we violate the equality (Q) and this 
can result in large magneto-asymmetry in the propaga- 
tion of unpolarized light. Such a large asymmetry, which 
we would refer to as magnetic field reversal asymmetry 
(BRAS in short), is induced by selectively applying the 
control field so as to break the time-reversal symmetry. 
It is important to note that we work with electric dipole 
transitions only. The control field is used to mimic the 
effects which occur in nature due to a combination of 
higher order multipole transitions. 

The structure of our paper is as follows. In Sec. II, 
we will discuss how one can use a control field to create 
large BRAS. We present a very simple physical model. 
We present relevant analytical results. In Sec. Ill, we 
introduce another model, where the control field is ap- 
plied such that we get a ladder system. We present the 
analytical results for BRAS in such a system. We also 
discuss the effect of atomic motion on BRAS in Sec. IV. 



II. LARGE MAGNETIC FIELD REVERSAL 
ASYMMETRY USING EIT 



The application of a coherent pump leads to the well- 
studied electromagnetically induced transparency (EIT) 
jl5| and coherent population trapping (CPT) pq . The 
usage of pump and probe in a Lambda configuration is es- 
pecially useful in suppressing the absorption of the probe, 
particularly if the lower levels of the A configuration are 
metastable and if the pump is applied between initially 
unoccupied levels. We explore how EIT can help in pro- 
ducing large BRAS. We first explain the basic idea in 
qualitative terms and then would produce detailed re- 
sults using density matrix equations for several systems 
of interest. 

Consider the following scenario. Let us consider first 
the case when B is applied parallel to the direction of 
propagation of the electromagnetic field. Suppose the 
control field is applied such that the <r+ component be- 
comes transparent, i.e., Im[x+(-B)] ~ 0. For magnetic 
field bigger than the typical linewidth,the component cr_ 
is off- resonant. Thus <r_ exhibits very little absorption 
Im(x-) ~ 0. Under such conditions, Eq. (||) shows that 
the transmitted intensity « /. Now if the direction of the 
magnetic field is reversed, then we easily find the situa- 
tion when <r_|_ component becomes off resonant from the 
corresponding transition, i.e., Im(x+) ~ 0; the tr_ com- 
ponent can become resonant and suffers large absorption, 
i.e., exhibits large Im(x-). This gives rise to an intensity 



~ 1/2. Thus the transmittivity reduces by a factor 1/2 
upon reversal of the magnetic field. It is thus clear how 
coherent fields can be used to create large BRAS. 

32p i/2 |B_ le > >= | 

F = 2 1 l i " A 

A 4 

' Q. 



3 2 S 



X 



=0 



1/2 
F=2 

lg> 



C 



-3B 



=+1 



lf> 



FIG. 1. The Na hyperfine level configuration is shown 
here. Here, B > is the applied magnetic field strength, 
2g± are the probe Rabi frequencies for the a± components, 
and Q is the half of the pump Rabi frequency. The respec- 
tive detunings 5 — [ui p — ui e+g (B — 0] and A for the probe 
and pump fields are defined with respect to the energy sepa- 
ration between the levels (|3 2 Pi/2;F = 2,m,F = 0), \g)) and 
(|3 2 P 1/2 ;F = 2,m F = 0), \3 2 S 1 / 2 ; F = 2, m F = 0}), respec- 
tively. Changing the direction of the magnetic field inter- 
changes the positions of |e_) and |e+). Besides the level |/} 
moves above the dashed line for |3 2 S , i/ 2 ; F — 2, ttif = 0). 

We now demonstrate the feasibility of these ideas. We 
consider a configuration [see Fig. [I]] which can be found, 
for example, in hyperfine levels of 23 Na |17|. The level 
\g) (\3 2 Si/ 2 ',F = l,mp = 0)) is coupled to the up- 
per level |e_) {\3 2 P 1/2 ;F = 2,m F = -1)) and |e+) 
(|3 2 Pi/2;^ = 2, nip = +1)) by the <7_ and a + com- 
ponents of the probe field. The susceptibilities for the 
two components of the probe acting on the transitions 
\g) <-* |e_) and \g) <-> |e+) are given by 



X-(B) 
X+(B) 



-i<ya 



i(S-B)-T e _ g 



i(S + B)-T. 



(9a) 

,5 = uj-Lj e+ g(B = 0), (9b) 



e+9 



where, «o is given by N\d\ 2 /Ti'y and is related to the 
absorption in the line-center for B = 0. It should be 
borne in mind that B represents the Zeeman splitting of 
the level mj? = —1. Thus B has the unit of frequency. 
Here 2j is the spontaneous decay rate from the level 
|e_), T e _ g = 7(r e+s = 47/3) is the decay rate of the 
off-diagonal density matrix elements between level |e_) 
(|e + )) and \g), N is the atomic number density, |d| is the 
dipole moment matrix element between the levels |e_) 
and \g), and 5 is the detuning of the probe field from the 
\g) \3 2 Pi/ 2 ',F — 2,m F = 0) transition. Note that 5 
would always be defined with respect to the levels in the 
absence of the magnetic field. Using the Eqs. (||) and (||), 



2 



one easily finds that the relation (ph holds for all 8. In all 
equations, B would be considered as positive quantity. 

To create a large asymmetry between the output inten- 
sities Io(B) and Io(—B), we now apply a coherent control 
field 



Ep(z,t)=E p (z)e- iu » t + < 



(10) 



on the transition \e + ) |/) (\3 2 Si/ 2 ', F = 2, rap = +1)). 
This modifies the susceptibility x+ of the a+ component 
to 



X+(B) 
A 



-i iao [i{6- A + 3B)-T fg ] 



[»(*• 



B)-T e+g ][i(8 
f(B = 0). 



a + w) -r 



fg\ 



M 2 ' 
(11) 



Here, A = 2B is the detuning of the pump field from 
the transition \3 2 P 1/2 ;F = 2,m F = 0) <-> |3 2 Si/ 2 ;.F = 

2,itif = 0) transition [see Fig. |J, fl — d e+ f.£ p /h is the 
half of the pump Rabi frequency The parameter T f g rep- 
resents the collisional dephasing between the states |/) 
and \g). In what follows we use T f g — 0. The level |/) is 
Zeeman separated from the level |3 Sim', F = 2, nip = 0) 
by an amount of 3B, whereas the levels \e±) are sepa- 
rated by an amount =F_B. These can be calculated from 
the Lande-g factor of the corresponding levels. The sus- 
ceptibility X- remains the same as in (J9|a) . Note that in 
presence of the control field, the response of the system is 
equivalent to a two-level system comprised of (|e_), \g)) 
[for er_ component] and a A-system [for er + component] 
comprised of (|e+), |/), \g)) connected via the common 
level \g). 
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FIG. 2. The variation of imaginary parts of the suscepti- 
bilities x+ (solid curve) and X- (dashed curve) in units of 
oto with probe detuning 5/7 are shown for the x±{B) [(a)] 
and X±(—B) [(b)]. The parameters used here are Q = O.57, 
B — 57 corresponding to 105 G, r e+9 = 47/3, Tf g = 0, 
r e _ 9 = 7, and A = 25. 

It is clear that applying a coherent pump field, one 
can generate an EIT window at 5 = — B (cf., A = 2B) 
for the ct+ component [Im (x+) — 0]. On the other 
hand, the absorption peak of the er_ component occurs 
at 5 = B. Thus, this component suffers a little absorp- 
tion [Im (x~) ~ 0] at 5 = —B as the field is far detuned 
from the |e_) <-» \g) transition as long as we choose the 
magnetic field much larger than the width of the tran- 
sition [see Fig. 0(a)] . Thus the unpolarized probe field 



travels through the medium almost unattenuated. The 
transmittivity T(B) = Iq(B)/I becomes almost unity at 
5 = — B as obvious from Eq. (|f) [see Fig. B], 




20.0 



FIG. 3. This figure shows the variation of the transmittiv- 
ity T(B) (solid curve) and T(—B) (dashed curve) with re- 
spect to probe field detuning 5/7. The parameters used here 
are N = 10 10 atoms cm -3 , A = 589 nm, and L = 1 cm. All 
the other parameters used are the same as in Fig. ^. 

If now the direction of the magnetic field is reversed 
(B — > —B), then the corresponding susceptibilities for 
<7_ and cr+ polarizations become 



X-(-B) 



i(5 



X+(-B) = 



B)-T e _ g 1 
-ija [i(5 - 



(12a) 



A — 3B) — r 



[i(8-B)-T e+g ][i(5-A-3B)-T 



(12b) 



We continue to take the quantization axis as defined by 
the direction of the propagation of the electromagnetic 
field. Clearly now at 8 = —B, X-(~B) has absorption 
peak and cr_ component of the probe will be absorbed. 
If we continue to use A = 2B, i.e., if we keep the con- 
trol laser frequencies fixed while we change the direction 
of the magnetic field, then x+{— B) exhibits resonances 
at 8 = 3B ± V '4B 2 + f2 2 , both of which are far away 
from the point 8 = —B unless we choose Vt 2 = 12B 2 . 
Clearly, for 8 = —B and tt 2 ^ 12B 2 , the a + component 
of the probe will suffer very little absorption. This is 
in contrast to the behavior of the cr_ component which 
will be attenuated by the medium. Thus the output field 
would essentially have the contribution from the cr+ com- 
ponent. The transmittivity T(-B) = I (-B)/I of the 
medium decreases to about 1/2. Thus by using EIT we 
can produce the result T(B) w 2T(—B), i.e., we can alter 
the transmittivity of the medium by just reversing direc- 
tion of the magnetic field. The equality (||) no longer is 
valid and the medium behaves like a chiral medium. This 
becomes quite clear from the Fig. |[ at 8 = —B. 
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A quite a different result is obtained by choosing the 
parameter region differently. For the choice of the exter- 
nal field strength fl = and A = 25, 5 = -B, 
the ct+ component gets absorbed significantly if the di- 
rection of the magnetic field is opposite to the direction 
of the propagation of the field. Thus T(—B) becomes 
insignificant compared to T(B) as shown in Fig. ^. For 
larger values of B, the result is shown in the inset. Here 
T{B)/T(-B) is in the range 2 to 3. Clearly, the case 
displayed in Fig. ^ is quite an unusual one. Such a large 
asymmetry in the dichroism of unpolarized light is the 
result of the application of a coherent control field whose 
parameters are chosen suitably. 



Im[5c + (-B)] 

5 x jmgjg)] 




FIG. 5. The variation of the imaginary parts of the suscep- 
tibilities x+{~B) and X-(-B) in units of Qo with the magnetic 
field B/"/ for 5 = — B. The inset shows the variation T(—B) 
with B in the vicinity of Q = 2y/3B. The parameters used 
here are the same as in Fig. ^j. 



III. LARGE MAGNETIC FIELD REVERSAL 
ASYMMETRY IN LADDER SYSTEM 



The behavior shown in the Fig. [I| is easily unterstood 
from the magnitudes of the imaginary parts of the sus- 
ceptibilties ^±(±-8). In the parameter domain under 
consideration, lm[x+(B)} = (EIT); Im[x-(-B) = a , 
because the ct_ component is on resonance for B an- 
tiparallel to the direction of propagation. Further as 
shown in the Fig. [|, in the region around ft = 2y/3B, 
Im[x-(-B)] -C Im[x+(— B)]. Thus both cr_ and a + com- 
ponents are absorbed if B is antiparallel k, making the 
medium opaque as shown in the inset of Fig. 5. The 
opacity disappears if direction of B is reversed (see Fig. 
4, solid curve). For values of B away from the equality 
il = 2V3B, Im[x+(— B-\] decreases leading to an increase 
in the transmission T(—B). 
0.5 




FIG. 4. The transmittivities T(B) and T(-B) calculated 
at the value 8 = —B is plotted here with respect to B/7, for 
f2 = IO7 and A = 2B. The inset shows the magnetic field de- 
pendence of the ratio T(B)/T(—B) for the same parameters. 
All the other parameters are the same as in Fig. pi 



It may be recalled that there are many different situ- 
ations where pump cannot be applied in a Lambda con- 
figuration. This, say, for example, is the case for 40 Ca. 
The relevant level configuration is shown in Fig. |^ [Q. 
The level \g) (|4s 2 ;j = Q,rrij = 0)) is coupled to |e+) 
(|4s4p;j = = +1)) and |e_) (|4s4p;j = l,rrij = 

— 1)) via the a+ and cr_ components of the input un- 
polarized probe field, respectively. In this configuration, 
the susceptibilities of the two circularly polarized com- 
ponents of the probe are given by 



X+[B) = 
X-(B) = 



i(5-B)-Y 
i(6 + B)-j' 



(13a) 

6 = uj -uj e _ g {B = 0), (13b) 



where, «o = N\d\ 2 /h^, N is the number density of the 
medium, 5 is the detuning of the probe field with respect 
to the \g) <-> |0) (|4s4p;j = l,TOj = 0)) transition, |rf| 
is the magnitude of the dipole moment matrix element 
between the levels |e+) and \g), and 27 is the decay rate 
from the levels |e + ) and |e_) to the level \g) . Note that 
the imaginary parts of the susceptibilities (O) are peaked 
at S = B and 5 = —B, respectively, and clearly predict 
perfect symmetry in the transmittivity of the medium 
upon reversal of the direction of the magnetic field. 
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FIG. 7. The variation of imaginary parts of the suscepti- 
bilities x+ (solid curve) and \- (dashed curve) in units of ao 
with respective to the probe detuning 5/7 are shown here for 
X±{B) [(a)] and X±(—B) [(b)]. The parameters used here are 

= O.57, B = 57 corresponding to 123 G, T = 0.457, and 
A = -B. 
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FIG. 6. The 40 Ca level configuration has been shown here. 
Here, 2g± are the probe Rabi frequencies for the a± compo- 
nents, B is the magnetic field strength, Q. is the half of the 
pump Rabi frequency, 8 and A are the respective detuning 
for the probe and the pump fields. These detunings are de- 
fined with respect to the energy separation between the levels 
(|0),| ff » and (| 0), I/)), respectively. 

We will now show how one can use a coherent control 
field to create asymmetry between T(B) and T{—B). We 
apply a coherent pump (CLO) to couple |e+) with a higher 



excited level |/) (\4p ;j 
quency 2f2 = 2df e+ .£ p /h. 
ponent now changes to p 



= 0, rrij — 0)) with Rabi fre- 
The susceptibility for <r+ com- 



X+(B) = 
A = 



-ija [i{A + S) 



[i(5-B)-~/][i(A + 8) 

J P ~ ^/e+C 8 = 0), 



(14) 



where, T — 0.5(X e+g /\f e+ ) 3 j = 0.457 is the spontaneous 
decay rate of the upper level |/) [cf., X e+g = 422.7 nm 
and A/ e+ = 551.3 nm], X a /3 is the wavelength of the tran- 
sition between la) and \/3), A = — B is the detuning of 
the pump field from the |/) «-> |0) transition [see Fig. |). 
Thus a transparency dip in the absorption profile of the 
er_|_ component at 5 = B is generated and the er_ compo- 
nent remains far-detuned from the corresponding transi- 
tion, as shown in Fig. 7|(a). Note that the transparency 
for cr + is not total which is in contrast to a Lambda 
system. We display in the Fig. || the behavior of the 
transmittivity T{B) of the medium as a function of the 
detuning. 



Now upon reversal of the magnetic field direction, the 
cr + component gets detuned from corresponding transi- 
tion and thereby suffers little absorption. The er_ compo- 
nent, being resonant with the corresponding transition, 
gets largely attenuated inside the medium. This is clear 
from the Fig. |?](b). Thus the contribution to the trans- 
mittivity T(—B) comes primarily from the o+ compo- 
nent. The Fig. || exhibits the behavior of T(—B) as the 
frequency of the probe is changed. In the region of EIT, 
T(B) is several times of T(-B). 
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FIG. 8. The variation for the transmittivities T(B) (solid 
curve) and T(—B) (dashed curve) with the probe detuning 
5/7 are shown in this figure. The parameters used here are 
TV = 10 10 atoms cm" 3 , A e+9 = 422.7 nm, L = 1 cm, and the 
other parameters are the same as in Fig. ^ 




In Fig. g, we have shown how the ratio T(B)/T(-B) 
calculated at 5 = B is modified with change in the mag- 
nitude of the applied magnetic field for different control 
field Rabi frequencies. Note that for large B and fi, this 
ratio approaches the value of two, though for intermedi- 



ate values it can exceed two. 
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FIG. 9. The variation of the ratio T(B)/T{-B) calculated 
at 8 — B with the magnetic field strengths B/j are shown in 
this figure for different values of Q. All the parameters are 
the same as in Fig. ^. 



IV. MAGNETIC FIELD REVERSAL 
ASYMMETRY IN THE PROPAGATION OF AN 
UNPOLARIZED BEAM THROUGH A DOPPLER 
BROADENED MEDIUM 

We next consider the effects of Doppler broadening on 
the BRAS in a Lambda configuration. We would like 
to find parameter regions where T(B) and T(—B) could 
differ significantly. We identify the spatial dependence of 
the pump £ p (z) = e lkpZ and A„ = 2B + k p v z , where v z 
is the component of the atomic velocity in the direction 
of the propagation of the electric fields. We assume that 
the pump field propagates in the same direction k p as the 
probe field wave vector k and we further take k and k p 
to be approximately equal. 

We calculate the Doppler-averaged susceptibilities 
through the following relation: 



(X±(v z )) 



X±{v z )(JD{vz)dv z 



(15) 



This leads to the following approximate results: 



(X±(v z )) = 



7"o 



k(v± — iujj) 



(18) 



where, 



i{[i(S + B) 



r e+9 ]p + |tt| 2 }. 



kp 



P = i{8 



3B) 



fg> 



(19) 



i[i{5-B)-T e _ g ] 



Here we have used the expressions (||a) and ( |lT|) for 
X±(v z ). These susceptibilities ([l8]) are used to calcu- 
late the transmittivities at the point 8 = —B. In Fig. [h] 
we have shown the corresponding variation of T{B) and 
T{-B) with 5/ j. We find that the ratio T(B)/T(-B) 
increases to a value ~ 1.6, for a 6 cm medium. 
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FIG. 10. The variation in T(B)/T(-B) with the probe de- 
tuning 5/7 in a 6 cm Doppler-broadened medium. All the 
other parameters are the same as in Fig. ^| Note that T(B) 
is the transmission for B parallel to k. 

It is clear that, if we choose a longer medium in this 
case, the BRAS will be further enhanced. We have actu- 
ally also carried out numerically the integration ( |l5| ) . For 
the parameters of the Fig. [l(], the results do not change 
substantially. 



where, 



od{v z ) 



-v 2 z /2u 



(16) 



is the Maxwell-Boltzmann velocity distribution at a tem- 
perature T with the width cub = yj KbT/M, Kb is the 
Boltzmann constant, and M is the mass of an atom. 

The integration Jig ) results in a complex error func- 
tion pH . However, to have a physical understanding, 
we integrate ( |l5| ) by approximating crjj by a Lorentzian 
ol{v z ) of the width Cob = 2wcln2 Jl9[: 



pl{v z ) 



lob/v 



(17) 



V. CONCLUSIONS 

In conclusion, we have shown how one can make use of 
a coherent field to create large BRAS in the propagation 
of an unpolarized light. We have discussed two different 
system configurations. We have shown how EIT can be 
used very successfully to produce large BRAS. We have 
also analyzed the effect of Doppler broadening and found 
the interesting region of parameters with large asymme- 
try. 
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